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A  new  three  dimensional  model  of  very  low  frequency  wave  propagation 
in  the  earth-ionosphere  waveguide,  in  the  presence  of  localized  and 
transient  disturbances,  has  been  developed.  The  model  is  generally 
applicable  to  disturbances  of  different  shape  and  ongin,  although  the 
application  up  to  now  has  been  to  interpret  observations  of  lonosphenc 
disturbances  associated  with  lightning-induced  electron  precipitation  frorP 
the  magnetosphere.  This  work  constituted  the  PhD  disssertation  of  Mr. 
William  Poulsen. 
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This  project  involved  the  development  of  a  new  three  dimensional  computer-based  model  of 
subionospheric  VLF  propagation  in  the  presence  of  localized  disturbances  of  the  D-region. 
A  list  of  publications,  talks,  and  PhD  thesis  developed  under  this  contract  are  provided  in 
the  attachment  In  this  report,  we  briefly  summarize  the  program  highlights  with  reference 
to  the  listed  papers. 

The  first  year  of  the  three  year  project  was  devoted  to  the  development  of  a  three 
dimensional  scattering  code  based  chi  a  formulation  Erst  published  by  Wait  [1962].  An 
important  milestone  in  the  first  year  of  the  project  (inititu^  on  1  October  1987)  was  the 
acquisition  of  the  existing  two  dimensional  pn^agation  code,  the  Long  Wave  Propagation 
Capability  (LWPC),  developed  and  used  by  the  Naval  Ocean  Systems  Center  (NOSC). 

This  NOSC/LWPC  code  was  acquired  in  July  1988  through  courtesy  of  Dr.  J.  Ferguson  of 
NOSC. 

An  initial  paper  describing  the  basic  aspects  of  the  new  three  dimensional  formulation  and 
its  application  to  special  cases  involving  single  dominant  waveguide  modes  was  submitted 
to  JGR  during  the  second  year  of  the  program  and  was  published  early  in  1990  [Poulsen  et 
al.,  1990a]. 

The  integration  of  the  three  dimensional  scattering  model  into  the  propagation  code  so  as  to 
allow  the  treatment  of  cases  with  multiple  dcxninant  modes  was  carri^  out  during  the  third 
year  of  the  program,  and  the  first  application  of  the  new  methodology  was  reported  in  the 
Proceedings  of  the  Ionospheric  Effects  Symposium  (lES)  [Poulsen  etal.,  19Wb]. 

A  fundamental  assuiiq)tion  adopted  in  the  develq>ment  of  the  new  three  dimensional  model 
was  to  neglect  coupling  between  wavegtude  modes  within  the  localized  scattering  center. 
This  assumption  was  shown  to  be  fully  justified  for  most  types  of  ionospheric  d^turbances 
that  are  existed  to  be  produced  by  li^tning-induced  electron  precipitation.  However,  due 
to  this  assuiziption,  the  applicability  of  the  new  model  is  limited  to  disturbances  with 
transverse  dimnsions  of  at  least  a  few  times  the  wavelength  of  die  scattered  waves.  A 
preliminary  calculation  of  the  scattering  from  very  small  (~S  km)  distinbances  indicated  that 
much  higher  electron  precepitation  fluxes  would  be  required  to  produce  observed 
subionospheric  VLF  amplitude^hase  changes  [Bell  et  al.,  1990].  The  conq)lete 
generali^cxi  of  the  model  so  that  it  can  be  applicable  for  the  treatment  of  ionospheric 
disturbances  of  arbitrary  shape  was  left  as  a  future  objective. 
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All  aspects  of  this  project  was  carried  as  part  of  (or  in  connection  with)  the  RiD 
Dissertation  of  Mr.  W.  L.  Poulsen.  Although  the  funding  for  the  project  terminated  on  30 
September  1990,  the  completion  of  the  PhD  Thesis  and  associated  research  papers  were 
unavoidably  flayed  beyond  that  time.  Mr.  Poulsen  has  now  successfully  defended  his 
PhD  thesis  on  7  Janu^  1991  (see  enclosed  Abstract)  and  is  currendy  completing  his 
written  dissertation  with  final  submittal  and  graduation  expected  in  April  1^1. 

At  least  two  research  papers  are  now  being  prepared  based  on  the  PhD  Thesis  for 
submission  to  RouUo  Science  and  Journal  of  Geophysical  Research.  The  titles  of  these 
papers  are  provided  below.  These  manuscripts  are  planned  to  be  completed  and  submitted 
to  the  journals  by  the  end  of  April  1991. 
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